Kolmogorov scaling from random force fields 
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We show that the classical Kolmogorov and Richardson scaling laws in fully developed turbu- 
lence are consistent with a random Gaussian force field. Numerical simulations of a shell model 
approximation to the Navier-Stokes equations suggest that the fluctuations in the force (acceler- 
ation) field are scale independent throughout the inertial regime. We conjecture that Lagrangian 
statistics of the relative velocity in a turbulent fiow is determined by the typical force field, whereas 
the multiscaling is associated to extreme events in the force field fiuctuations. 



In studies of fully developed turbulence, two discov- 
eries are highly noticeable as fundamental and seminal. 
One regards L.R. Richardson's study of the enhanced dis- 
persion of particles advcctcd by a turbulent flow [l| . The 
other result is Kolmogorov's fundamental derivation, es- 
sentially based on dimensional arguments, of the energy 
spectrum in fully developed turbulence . Both theories 
employ the energy cascade, from the integral scale down 
to the dissipation scale, as the paradigmatic physical pic- 
ture of the energy dissipation flow. Indeed, pair-particles 
passively advected by turbulence exhibit a superdiffu- 
sive behavior with their relative distance given by the 
Richardson's scaling Q , [3| . In contrast to the pair dis- 
persion, the motion of a single particle is determined by 
the correlation time of the underlying velocity field, such 
that it is transported ballistically for times smaller than 
the correlation time and diffuses normally for larger time 
scales m. In the cascade scenario, the pair-particle su- 
pcrdiffusion is due to the large energy jumps between the 
eddies containing each particle. However, in the velocity 
space no superdiffusive behavior is needed to substanti- 
ate this jumpy motion. 

In this Letter, we show that the velocity increments 
generated by a white-noise force field are sufficient to 
generate the superdiffusive behavior, as well as the Kol- 
mogorov energy spectrum. To put it in very simple 
terms: integrating 'up' from the random acceleration 
field to the velocity field and subsequently to the dis- 
placement is enough to reproduce the well-known scaling 
laws. 

To clarify the underlying physical picture, we consider 
a simple stochastic model of relative dispersion in a white 
noise acceleration field given by 
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where Av{t) = vi{t) ~ V2{t) is the velocity difference be- 
tween the two particles moving along the two trajectories 
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FIG. 1: Two particles being advected in a random force field, 
generated by the GOY shell model. 



ri{t) and r2{t), and AF{t) = Fi{t) - F2{t) is the relative 
force. The prcfactor 4 instead of the usual factor 2 ap- 
pearing in the force correlation is due to the parametriza- 
tion of the relative dispersion in terms of the diffusion 
constant e* for a single particle in the velocity space. 
The 5— function may have, in principle, a finite width 
given by the time correlation of the relative random field 
along the two trajectories. This width will be determined 
both by the time it takes to pass a correlation length for 
a given force realization, and the time it takes to change 
the force in a certain point of the system. 

In this set up, the relative velocity field, Av{t) = 
J* AF{s)ds, is a Wiener process with a Gaussian dis- 
tribution, namely 



{Av{ti)Av{t2)) =4e* min{ti,t2), 



(3) 



while the relative separation is described by a non- 
Gaussian distribution with the second moment satisfying 
the Richardson's scaling, that is 



4e* 



-t' 



(4) 
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FIG. 2: The squared relative acceleration AFAF and its in- 
finite moment versus Ar^. The thin lines are for the La- 
grangian trajectories where distances and accelerations are 
parameterized by the time of advection. The squares repre- 
sent the corresponding Eulerian measures of the same quan- 
tities. 

ity and distance, wc obtain the exact Kolmogorov scaling, 
{Av^it)) = 48i/3e*2/3(Ar2(i))i/3^ (5) 

in the Lagrangian framework (for the higher moments 
see @). Thus, Kolmogorov scaling is consistent with the 
assumption that the dispersion is driven by sufficiently 
random and uncorrelated acceleration fields. 

In deriving Eq. ([5]) we assumed that the relative veloc- 
ity is obtained by following the Lagrangian trajectories, 
which in a real turbulent flow may differ from the typical 
velocity increments separated by the distance r (Eulerian 
measurement of the velocity differences) Q- 

Eq. m implies that 2e* is the diffusion constant for the 
relative velocity. For a Lagrangian stochastic flow gen- 
erated by the white noise acceleration field, e* can be 
estimated as 

{Av{t)AF{t)) = {[ dsAF{s) ■ AF{t)) 
Jo 

= I ds{AF{s)- AF{t)) 
Jo 

^ 12e [ S{t~s)ds^ I2e*, (6) 
Jo 

where the additional factor 3 compared to eq.[3]is related 
to the 3D system. From dimensional considerations, e* 
has the same units [length'^ /time^] as the standard en- 
ergy dissipation rate e characterizing the turbulence cas- 
cade. 

To examine how the Lagrangian white-noise accelera- 
tion relates to the anomalous scaling laws in a more re- 
alistic turbulent field, we consider the kinematics of pair 
particles advected by the homogeneous turbulent flow ob- 
tained by a real-space transformation of the GOY shell 



model [1, This model proposed originally by Glcn- 
zer, Yamada and Ohkitani [lO, lllll provides a description 
of the turbulent motion embodied in the Navier-Stokes 
equations. The GOY model is formulated on a iV-discrete 
set of wavenumbers, fc„ = 2", with the associated Fourier 
modes m„ evolving according to 

^K-lK-2) + fSn,U (7) 

for 71 = 1 • • • A^. The coefficients of the non-linear 
terms are constrained by two conservation laws, namely 
the total energy, E = IwnP, and the helicity (for 
3d), H = J2n{—^)^kn\un\, 01 the enstrophy (for 2d) , 
Z — A:^|u„p, in the inviscid limit, i.e. f = v = 
Therefore, they may be expressed in terms of a free pa- 
rameter only i5 S [0,2], a„ = I, bn+i = —S, Cn+2 = 
— (1 — (5). As observed by Kadanoff [13], one obtains 
the canonical value e = 1/2, when the 3d- helicity is 
conserved. The set ([7]) of A^-coupled ordinary differen- 
tial equations can be numerically integrated by standard 
techniques [3]. We have used standard parameters in 
this paper TV = 19, = f 0"^, fco = 2 • lO'^, / = 5 • IQ-^. 

The GOY model is defined in A:-space but we study 
particle dispersion in direct space obtained by an inverse 
Fourier transform @ of the form 

N 

v{r, t)^Y.'^^ (^)e''" + c- c] . (8) 

n=l 

Here the wavevectors are fc„ = knCn where Cn is a unit 
vector in a random direction, for each shell n and c„ are 
unit vectors in random directions. We ensure that the 
velocity field is incompressible, V • = 0, by constraining 
Cn ■ e „ = 0, Vn. In our numerical computations we con- 
sider the vectors c„ and e„ quenched in time but averaged 
over many different realizations of these. 

As an example of the motion in this field. Fig. [T] shows 
the trajectories of two passively advected particles. As 
the relative distance diverges in time, the two particles 
experience different force fields, which in turn typically 
increase the difference in the relative velocities of the 
two particles. The figure shows the individual particles 
as they are advected, first together and later diverging 
away from each other when they are encased in different 
eddies. 

Fig. [2] examines the noise in the effective force field 
(AFAF) for the relative motion of the two advected par- 
ticles. In Fig [2^) wc use viscosity v = I0~^, with a Kol- 
mogorov scale Ar ^ I.O • 10"'*. The noise amplitude is 
plotted versus the average square distance between the 
particles (Ar^) — {{vi{t) — r2(t))^) , with the time as 
parametrization of the curves, as in eq. [S] The average 
is over independent trials of the two advected particles. 
One observes that both the typical noise and the max- 
imum value at any distance is constant throughout the 
inertial range, i.e. above the Kolmogorov scale. 
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FIG. 3: The squared relative velocity (Au)^ and its infinite 
moment versus (Ar)^ (ie a representation of the structure 
function). The thin lines are for the Lagrangian trajectories 
where distances and velocities are parametriced by the time 
of advection. The filled squares represent are the correspond- 
ing Eulerian measures. The straight line represents standard 
Kolmogorov scaling {Av'^) oc (Ar^)^^^. 



We conclude that the force field is equivalent to Gaus- 
sian white noise, and therefore the structure function of 
this turbulent field should be close to the one predicted 
by Eq. ([5]). This is confirmed in Fig. [3]where we show the 
deviations in velocity as a function of the square distance 
between the particles. One indeed sees that {Av^{t)) 
versus {Ar^{t)) scales with an exponent close to 1/3 in 
agreement with our expectations. For completeness, we 
also show the infinite moment of the velocity, and remark 
that this higher moment scales with an exponent close to 
0.23. This signals multidiffusion 15[ where extreme ve- 
locity differences sometimes, but rarely, are reached after 
short separations. In the current context, we see these 
extreme deviations as a measure of very unlikely and in- 
termittent events which only add little to the typical be- 
havior of the flow. Indeed also the Eulerian statistics 
shows clear multiscaling as expected Q. 

While our intuition has been based on the Lagrangian 
picture of advected particles, it is remarkable that 
the corresponding Eulerian quantities behaves similarly. 
This is demonstrated in simulations where we now fix the 
distance between two points, and then calculate respec- 
tively the difference in velocity and acceleration. The 
continuous curves in Figs. [2] and [3] show how (Ai^^(r)) 
and (Aw^(r)) vary with the square relative distance be- 
tween the investigated points. From Fig. [2] we see that 
the value of the plateau for the random force field is a di- 
rect consequence of its random expectation at any large 
distance. Therefore, there is nothing special about the 
selection of advected points in the Lagrangian case. In 
fact the onset of the plateau is slightly sharper in the Eu- 
lerian case, presumably refiecting averaging associated to 
the underlying time parameter in the Lagrangian advec- 
tion. Similarly, there is no significant difference for the 
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FIG. 4: A measure of an effective difi^usion constant in ve- 
locity: {AvAF) versus Ar^. As in the previous figures the 
thin line describes Lagrangian trajectories parametriced by 
the time of advection. The filled squares correspond to the 
Eulerian case. 



structure functions shown in Fig. [3] 

Using Eq. [H] we estimate {AvAF) ~ 0.1 through- 
out the inertial range in the GOY model simulations, 
see Fig. 01 This value of the effective velocity diffusion 
constant is larger than the average energy dissipation at 
the Kolmogorov scale, estimated from the energy input 
Re{ul ■ f) = 0.001 in the GOY model. This discrepancy 
in the effective diffusion terms we attribute to the huge 
contributions from the spikes in the acceleration which 
are absent in the simple white noise calculation of Eq. (5] 
These spikes also gives rise to multidiffusion, as discussed 
above. 

We believe that the acceleration field as shown in Fig. 
[5] should be experimentally accessible either by particle 
tracking in a 3-D flow [l^ or from probe measurements 
in channel flows employing the Taylor hypothesis. In 
the first case the acceleration is easily estimated from 
the temporal variations in the velocity field of the 3-D 
advected particles. 

Overall we have seen that the force field reaches an av- 
erage value that is independent on the distance between 
the advected points in the turbulent fluid. Already at 
distances slightly above the Kolmogorov scale the two 
particles often receive random "kicks" which are as large 
at small scales as they are at the integral scale. Thus, 
huge accelerations are associated to the very small scales, 
presumably to the core of eddies at the verge of their de- 
struction by dissipation. The acceleration between two 
particles are primarily dependent on how close each of 
them are to the center of an eddy. When examining the 
distribution of the accelerations at a fixed distance we 
observe a broad power law like behavior with a cutoff 
which is independent of the distance (as demonstrated 
by the constant max norm). The size of the cutoff is 
determined by the size of the forcing and the scale at 
which this forcing is acting (in our simulations, the scale 
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is Ar = 1). 

In conclusion, the motion associated to the relatively 
slow turn over dynamics of the large eddies is not needed 
for obtaining Richardson or Kolmogorov statistics. 
These two seminal laws are primarily a consequence of 
the random force field that fluctuates with an amplitude 
set by the system size and with a correlation time set by 



the Kolmogorov scale. 
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